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SECTION  I 

1.  0  TECHNICAL  REPORT  SUMMARY 

The  feasibility  of  using  seismic  waves  to  detect  discontinuities  in  rocks 
is  well  established.  The  capability  to  detect  and  map  such  discontinuities, 
whether  from  faulting,  changes  in  rock  type,  voids,  or  abandoned  mine 
workings,  would  be  useful  in  deep  underground  excavation  and  mining 
operations. 

At  the  present  time,  this  capability  exists  for  distances  of  only  a  few  feet. 

The  reasons  are  twofold: 

(1)  The  ultrasonic  frequencies  used  are  absorbed  too  rapidly 
in  the  rock,  and 

(2)  Unambiguous  interpretation  is  made  difficult  by  the  presence 
of  multiple  wave  arrivals. 

The  solution  to  the  first  problem  is  a  more  suitable  seismic  source.  We 
have  investigated,  analytically  and  experimentally,  the  characteristics  of 
several  possible  sources  for  use  on  the  tunnel  working  face  and  wall 
surfaces  and  also  in  small  diameter  drill  holes. 

The  lack  of  a  suitable  down-hole  source  led  to  the  selection  and  development 
of  an  underwater  spark  source  based  on  the  electrohydraulic  effect.  Although 
this  principle  had  been  previously  applied  as  a  source  for  ocean  acoustics  and 
offshore  seismic  exploration,  the  electrode  configurations,  energy  levels,  and 
seismic  pulse  characteristics  were  not  suitable  for  a  portable,  down-hole 
source.  The  device  developed  provides  an  estimated  12  joules  acoustic 
output  in  a  pulse  of  30-100  microseconds  duration. 
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Our  approach  to  the  problem  of  multiple  arrivals  is  the  use  of  seismic 
array  techniques  to  separate  the  reflections  and  unwanted  arrivals.  Two 
methods,  developed  during  a  previous  program,  are  to  be  verified  and 
evaluated  using  field  recorded  data  during  the  current  program.  These 
techniques  are: 

(1)  Beamforming,  which  uses  array  signal  processing  to 
enhance  arrivals  from  specific  directions,  and 

(2)  A  three-dimensional  CRT  display  where  displayed  intensity 
is  proportional  to  the  received  seismic  signal  at  a  point  in 
time  and  display  spatial  coordinates  are  proportional  to 
array  element  spatial  coordinates. 

Both  methods  utilize  crossed  array  configurations  which  can  be  conveniently 
implemented  in  confined  underground  spaces  with  down-hole  seismic  source 
and  receiver  arrays. 

To  provide  field  data  for  evaluation  of  the  array  methods,  a  portable  prototype 
seismic /acoustic  system  Has  been  developed.  Conventional  seismic  systems 
generally  record  on  magnetic  tape  the  signals  received  simultaneously  on  a 
large  array  of  receivers  for  subsequent  processing  and  display.  In  contrast, 
our  system  approach  is  to  present  the  seismic  signals  as  a  variable  intensity 
CRT  display  and  record  directly  as  density  variations  on  Polaroid  film.  The 
film  can  be  interpreted  immediately  in  the  field,  or  digitized  for  further 
processing. 

The  system  includes  a  pneumatically  clamped  borehole  receiver  and  the 
borehole  spark  source  referred  to  above.  Provision  was  also  made  for  an 
optional  pulsed  piezoelectric  transducer.  Nonrepetitive  sources,  such  as 
impacting  hammers,  can  also  be  used  with  the  system. 


In  operation,  the  single  repetitive  source  and  single  receiver  provide  one 
array  output  signal  at  a  time.  Each  signal  is  recorded  on  film  and  the 
receiver  repositioned  and  the  film  advanced  for  the  next  record.  In  this 
way,  an  array  consisting  of  as  many  as  512  elements  can  be  recorded  on 
a  single  4x5  Polaroid  film. 
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2.0  INTRODUCTION 

The  objectives  of  this  program  are  to  verify  and  evaluate  two  promising 
seismic  array  techniques  for  use  in  undergrr—  1  hard-rock  excav’  ins. 
The  two  array  techniques  to  be  evaluated 

(1)  A  three-dimensional  method  of  displaying  the  output 
signals  from  a  two-dimensional  array  of  transmitting 
sources  and  receivers  to  enable  an  operator  to 
visually  detect  reflections  from  discontinuities  within 
the  probed  rock  volume;  and 

(2)  A  beamforming  method  developed  in  a  previous  study 
which  employs  sophisticated  array  signal  processing 
to  enhance  reflections  returned  from  specific 
directions. 

The  evaluation  and  comparison  of  these  two  techniques  for  detection  and 
location  of  rock  discontinuities  are  directed  tow  ard  future  applications  of 
the  ARPA  Rock  Mechanics  and  Rapid  Excavation  Program.  Thus,  we 
are  primarily  concerned  with  deep  underground  openings  in  crystalline 
rocks.  Specifically,  we  would  antic  .  ate  relatively  few  reflecting 
discontinuities,  but  those  that  did  exist  would  be  of  primary  interest: 
changes  in  rock  type,  altered  zones,  or  faults.  Since  these  features  are 
also  of  interest  for  mining  operations,  successful  development  of  an 
underground  seismic  mapping  technique  might  also  be  expected  to  be 
useful  in  under  ground  mining  operations. 


The  underground  use  of  seismic  methods  is  viewed  as  complementary 
to  exploratory  drilling,  rather  than  as  a  replacement.  Although 
incapable  of  the  accuracy  provided  by  drilling,  seismic  methods  can 
economically  probe  substantially  larger  volumes  of  rock.  Discontinuities 
can  be  seismically  mapped  and  the  nature  of  the  discontinuity  determined 
from  a  single  borehole,  where  multiple  boreholes  might  otherwise  be 
required. 

Of  critical  importance  to  the  rapid  excavation  application  is  the  timeliness 
of  the  information  provided  by  the  seismic  system.  The  value  of  the 
information  is  determined  by  the  amount  of  advanced  warning  of  changing 
geologic  conditions  ahead  of  the  excavation  working  face.  This  requirement 
implies  that  the  seismic  system  approach  must  be  capable  of  being 
implemented  in  nearly  real  time  and  that  the  penetration  depth  be  as  large 
as  possible. 

The  attenuation  factor  for  seismic  waves  in  rock  is  proportional  to 
frequency.  Thus,  as  greater  penetration  is  required,  lower  frequencies 
must  be  used.  However,  since  the  ability  to  resolve  range  differences 
depends  on  the  wavelengths  employed,  there  will  be  some  transmitted 
pulse  spectrum  which  provides  the  best  compromise  between  resolution 
and  penetration. 

A  significant  part  of  the  present  program  is  concerned  with  seismic 
source  development  and  evaluation.  In  addition  to  the  frequency  content 
of  the  seismic  pulse,  other  source  considerations  include  pulse  duration, 
energy,  repeatability,  method  of  coupling  to  the  rock,  and  portability. 

A  range  of  at  least  50  meters  is  desirable  for  use  in  rapid  excavation 
applications.  Because  of  frequency-dependent  attenuation,  a  significant 
fraction  of  the  transmitted  pulse  energy  should  be  contained  below  5-10  KHz. 
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Moreover,  to  achieve  a  range  resolution  of  one  meter  for  a  typical  hard- 
rock  seismic  velocity  c*  5000  meters /second,  the  pulse  duration  should 
be  no  more  than  one  or  two  hundred  microseconds. 

The  seismic  pulse  can  be  applied  either  to  the  surface  of  the  rock  or  to 
the  interior  of  a  borehole.  Problems  arising  with  a  surface  source 
include. 

(1)  That  of  coupling  to  a  rough  and  uneven  surface, 

(2)  The  possibility  of  near-surface  fractures  caused  by 
blasting  or  rock  stress  release,  which  would  interfere 
with  propagation  of  the  seismic  pulse,  and 

(3)  Generation  of  interfering  Rayleigh  surface  waves. 

The  use  of  a  borehole  source  would  alleviate  these  problems,  but  would 
require  the  added  expense  of  drilling  the  hole  and  may  introduce  further 
complications. 

At  the  start  of  the  current  phase  of  the  program,  there  were  a  number  of 
surface*  type  sources  available  from  which  to  choose,  including  hammers 
and  other  impact  sources  as  well  as  piezoelectric  transducers.  Conse¬ 
quently,  effort  was  directed  toward  developing  a  suitable  down-hole  source 
and  evaluation  of  existing  surface  sources.  An  additional  requirement  was 
that  the  hole  diameter  required  for  the  source  be  no  more  than  about  1-1/2 
inches,  to  minimize  drilling  cost. 

Similar  considerations  apply  to  the  use  of  borehole  receivers.  The  same 
problems  arise  in  the  use  of  surface -emplaced  receivers  as  for  surface 
sources.  Hydrophones  or  pressure  sensors  in  a  water-filled  hole  are  an 
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alternative  possibility.,  but  the  presence  of  the  fluid  in  the  hole  allows 
other  wave  modes  to  propagate  along  the  hole.  These  waves  would 
tend  to  interfere  with  or  obscure  the  desired  refle^^ons.  Thus,  a 
second  goal  for  this  program  phase  has  been  the  development  of  a 
receiving  sensor  which  can  be  readily  clamped  to  the  wall  of  a  small 
diameter  borehole. 

A  third  goal  for  this  phase  has  been  the  development  of  an  optical 
recorder  which  displays  the  received  signals  as  a  CRT  intensity- 
modulation  and  provides  a  variable  density  film  record  of  the  seismic 
array  output.  The  three -dimensional  display  can  be  viewed  directly 
from  the  variable  density  negative,  which  can  also  be  digitized  for 
beamforming  purposes.  Usually,  in  seismic  recording  systems, 
magnetic  tape  is  used  as  an  intermediate  storage  medium.  However, 
for  simplicity  and  minimum  cost,  and  because  of  the  eventual  need  for 
real-time  viewing,  the  intensity -modulated  CRT  with  photographic  film 
recording  was  chosen. 

Together,  the  seismic  source,  receiver,  and  CRT  display -recorder 
constitute  a  relatively  inexpensive,  portable  seismic  array  system  for 
use  in  underground  excavations.  The  addition  of  an  optical  viewer  will 
allow  near  real-time  field  interpretation  using  the  three-dimensional 
display  method,  and  a  computer  module  for  beamforming  can  also  be 
added,  if  desirable,  at  some  future  date. 


3-1 


SECTION  III 

3.  0  SEISMIC  ARRAY  APPROACHES 

The  use  of  seismic  source/receiver  arrays  offers  several  advantages 
over  a  single  source  and  receiver  combination  for  reflection  seismology: 

•  Identification  of  particular  wave  types  is  enhanced. 

•  Reflected  waves  arriving  from  some  Darticular  direction 
can  be  enhanced  relative  to  other  interfering  waves. 

•  Direction  to  a  reflector  can  be  determined. 

Two  distinctly  different  array  concepts  will  be  evaluated  in  an  underground 
environment  during  the  next  phase  of  the  program.  These  two  approaches 
are: 

(1)  Beamforming  or  spatial  filtering,  and 

(2)  The  "three-dimensional  seismic  method". 

The  application  of  beamforming  to  the  hard  rock  excavation  problem  was 
investigated  under  a  previous  contract  (Ref.  7),  while  the  3-D  method  is 
a  new  way  to  present  seismic  array  data  which  has  been  recently  developed 
by  ESSO  Production  Research  Company  for  conventional  seismic 
exploration  (Refs.  3, 8). 


Regardless  of  the  array  concept  used  for  processing  and  display  of  the 

data,  a  two-dimensional  or  planar  array  is  required  by  the  three- 

dimensional  nature  of  the  underground  problem,  if  dip  and  strike  angles 

of  reflecting  surface  are  desired.  The  ideal  configuration  would  consist 

2 

of  a  large  square  array  of  N  seismometers  distributed  on  a  flat  surface 
(Fig.  3-1).  This  approach  would,  in  general,  be  impractical  for 
underground  use  because: 

(1)  A  sufficiently  large  flat  surface  would  probably  not  be 
available,  and 

o 

(2)  The  cost  of  providing  N  receiving  elements  would  be 
prohibitive. 

Alternatively,  a  single  scanning  receiver  could  be  used  with  a  repetitive 
source  (Fig.  3-1),  but  the  first  objection  would  still  be  valid.  The 
flatness  requirement  could  be  relaxed  somewhat  by  measuring  the 
displacement  of  each  receiver  position  relative  to  a  datum  plane  and 
providing  appropriate  time  delays  for  each  receiver  output  signal. 
However,  this  approach  also  seems  impractical  and  costly. 

An  alternative  is  to  use  crossed  linear  arrays.  Milder  and  Wells  (Ref.  5) 
have  described  a  crossed  array  consisting  of  a  line  array  of  receivers  and 
a  line  array  of  transmitters  for  acoustical  holography  (Fig.  3-2a).  With 
this  configuration,  only  2N  transducers  (N  receivers  and  N  transmitters) 
are  required  for  the  equivalent  of  a  square  array  of  N  receivers  [the 
(n,  m)^1  datum  results  from  listening  to  the  m^1  transmitter  with  the  n^1 
receiver]. 
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a.  Crossed  Array,  Stationary  Transmitter  b.  Crossed  Array, 

Array  Scanning  Transmitter 


Figure  3-2.  Alternatives  to  Square  Array 
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This  approach  was  also  adopted  for  the  3-D  method  for  seismic 
exploration  (Ref.  8  ).  There  are  several  advantages  of  the  crossed 
array  technique  for  underground  application.  For  example, 
perpendicular  holes  can  be  drilled  in  the  rock  for  the  two  line  arrays 
(Fig.  3-3).  The  inner  surfaces  of  the  holes  will  provide  reasonably 
smooth,  uniform  surfaces  for  coupling  the  seismic  energy  into  and  out 
of  the  rock.  The  holes  can  be  extended  to  any  practical  length  D,  so 
the  size  of  the  array  (hence,  angular  resolution)  is  not  restricted  by 
the  dimensions  of  the  mine  or  tunnel  opening.  The  holes  can  be 
drilled  into  the  walls,  floor,  and/or  roof  parallel  to  the  working  face, 
but  away  from  the  face  to  minimize  interference  with  excavation 
operations. 

The  crossed -array  data  can  also  be  synthesized  by: 

(1)  Using  a  single  scanning  transmitter  with  a  stationary 
receiving  array  (Fig.  3 -2b), 

(2)  Using  a  single  scanning  receiver  with  a  stationary 
transmitter  array,  or 

(3)  Using  a  scanning  receiver  and  a  scanning  transmitter. 

This  last  approach  for  data  generation  has  the  advantages  of  simplicity, 
low  cost,  and  ease  of  implementation,  and  has  been  selected  for  the 
evaluation  phase  of  the  program.  It  is  also  the  slowest,  so  that  we  might 
anticipate  one  of  the  other  approaches  being  more  suitable  for  implemen¬ 
tation  in  an  operational  system. 


A  crossed-array  configuration  was  also  recommended  for  beamforming 
as  a  result  of  the  previous  study  (Ref.  7).  The  configuration  investigated 
there  consisted  of  perpendicular  line  arrays  of  receivers,  with  a  single 
transmitter  at  the  intersection  of  the  two  lines.  This  configuration  could 
be  implemented  with  either  stationary  receiver  arrays  (Fig.  3 -4a),  or 
with  a  single  scanning  receiver  and  repetitive  source.  The  outputs  of 
each  array  were  separately  combined  to  form  beams  with  good  angular 
resolution  in  the  two  spatial  coordinates  normal  to  the  array  axis.  The 
beamformer  outputs  for  each  line  array  were  then  cross -correlated  with 
those  of  the  other  (Fig.  3 -4b)  to  provide  narrow  beams  in  three  spatial 
coordinates. 


3.2  Beamforming  Algorithms 

The  beamfor.ning  algorithms  previously  recommended  for  further  field - 
scale  investigation  were  as  follows: 

•  Digital  Multibeam  Steering  (DIMUS) 

«  DIMUS/Adaptive  Null 

These  recommendations  were  based  on  simplicity  of  implementation  and 
relative  performance  with  respect  to  other  possibilities  which  included: 

(1)  Delay -and-sum  beamforming 

(2)  Weighted  delay -and -sum  beamforming 

(3)  Fan  filtering 

(4)  Least-squares  (Wiener)  spatial  filtering 

(5)  Adaptive  beamforming 
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b.  Crossed  Correlation 
Figure  3-4.  Crossed  Array  for  Beamforming 
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The  DIMUS  and  DIMUS /Adaptive  Null  methods  both  employ  hard  clipping 
of  the  receiver  outputs.  That  is,  only  the  polarities  of  the  signals  are 
processed.  Clipping  provides  several  benefits  and  some  adverse  effects 
as  well.  Benefits  include: 

•  Simpler  implementation 

•  Insensitivity  to  amplitude  variations 

•  Strong  interferences  are  reduced,  relative  to  small  signals 

Simpler  implementation  results  from  processing  only  the  polarity  bit. 

This  reduces  analog-to-digital  conversion,  data  storage  and  digital 
processing  requirements.  Amplitude  variations  due  to  differences  in 
source /receiver  coupling  to  the  rock  across  the  array  or  variations  in 
source  output  are  unimportant  and  there  is  nc  need  for  AGC.  On  the  other 
hand,  disadvantages  include  the  broadening  of  narrow  pulses  with  tails 
(decrease  of  resolution),  and  the  possibility  cf  reducing  strong  signals 
relative  to  weak  interferences.  On  balance,  however,  the  advantages  of 
hard  clipping  appear  to  outweigh  the  disadvantages. 

3. 2. 1  The  DIMUS  Algorithm  —  DIMUS  beamforming  consists  of  delay  - 
and-sum  beamforming  on  hard -clipped  data: 

b(n)  =  E  sgn  [sk(n-mk)J 

where 

( 

j  +1,  z  a  0 

sgn(z)  =  ( 

-1,  z  <  0 
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til 

and  s^(n)  is  the  amplitude  of  the  signal  sensed  at  the  k  receiver  at 
the  n^  sample  instant  after  transmission  of  a  pulse  by  the  source.  The 
integer  mk  is  the  delay  applied  to  the  k^1  receiver  output  to: 

(1)  Steer  the  beam  in  a  particular  direction, 

(2)  Correct  for  wavefront  curvature  (near-field  focusing),  and 

(3)  Correct  for  array  out -of -planeness,  if  necessary. 

Fig.  3-5  is  a  block  diagram  for  a  generalization  of  DIMUS  b_dmforming 
in  which  weights  are  applied  to  each  receiver  output  to  modify  the 
beam  pattern. 

3. 2. 2  DIMUS  /Adaptive  Null  --  This  algorithm  allows  nulling  of  a  coherent 
off-axis  interference  which  arrives  coincidentally  with  an  on -axis  signal, 
using  hard-clipped  data.  The  interference  is  nulled  by  subtracting  from 
the  k**1  receiver  output  a  coherent  estimate  of  the  interference  obtained 
from  the  (k+l)st  receiver  output.  The  beam  output  is  given  by 

bta)  =  £  t f  k  (n-n'k)  '  wpk  ck+l  <n'mk'“pk>,;  PNo<n<  <P+1)N0-  1 
where 

ek(n)  =  sgn  [sk(n)J. 

The  integer  u  ^k  is  the  sample  delay  corresponding  to  the  propagation 
time  of  the  interference  between  adjacent  receiver  locations.  The  optimum 
value  of  jipk  is  chosen  by  maximizing 
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spk 


max 

i 


(p+l)N0-l 


Z 

n=pN0 


h  <n-mk>  ck+i  ‘“‘"V0 


It. 

on  the  p  sample  interval  of  length  NQ,  scanned  over  a  range  of  delays 
encompassing  the  maximum  interelement  propagation  time  for  interferences. 


The  weights  W  .  are  given  by 


W 


pk 


0*  1  MpjJ  £  6pk 

1;  I  Mpk'  >  6pk 


The  constraint  Wpk  =  0  if  |  n^\  i  6pk  prevents  nulling  of  a  signal 
within  the  beam,  provided  £pk  is  chosen  equivalent  to  a  half  beamwidth 
in  interelement  sample  delay. 


Figure  3-6  is  a  block  diagram  of  the  DIMUS /Adaptive  Null  beamformer. 


3.2.3  Cross -Correlation  (Time -Averaged  Product)  of  Beams  -Denoting 
the  two  crossed  arrays  by  their  axis  designation  (x,  y),  if  b.  (n)  is  the 
beamformed  output  of  the  iul  array,  the  cross  correlation  of  two  beams, 
which  provides  a  narrower  beam  than  either  of  the  individual  beams,  is 
given  by 


C(n) 


b  (n  + 1 )  b  (n  +  £). 
A  y 


The  correlation  window,  (2L+  1)  samples,  is  chosen  equivalent  to  a  pulse 
length. 

This  process  is  shown  schematically  in  Figure  3-4b. 
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3.  3  Three-Dimensional  (3-D)  Seismic  Method 

Three-dimensional  treatment  of  seismic  data  is  essential  for  underground 
applications.  In  the  beamforming  methods  described  above,  this  is 
accomplished  by  calculating  the  steering  delays,  m^,  to  correspond  to 
spherical  angle  pairs  (0,  <j>)  in  the  directions  of  interest.  After  applying  the 
appropriate  delays  to  all  receiver  outputs  and  summing,  coherent  arrivals 
will  add  with  similar  polarities  (i.e.,  in  phase).  The  spatial  phase 
coherence  of  reflected  waveforms  may  be  at  least  partially  destroyed  (i.  e. , 
wavefront  distorted),  by  one  of  several  causes,  including: 

•  Rock  inhomogeneity 

•  Nonplanar  reflecting  surface 

•  Nonplanar  array  locations 

•  Faulting 

As  a  result,  the  beamformer  output  will  be  reduced,  decreasing  the  likeli¬ 
hood  that  the  reflection  will  be  detected. 

The  three -dimensional  seismic  method  is  an  alternative  approach  which  is 
expected  +o  be  more  robust  for  the  conditions  listed  above.  The  3-D  method 
will  be  easier  to  implement,  but  will  place  a  greater  interpretation  burden 
on  the  operator. 

The  3-D  method  can  best  be  described  with  reference  to  a  single  source  at 
the  center  of  a  large  square  array  of  receivers.  The  outputs  of  the 
receivers,  in  response  to  a  transient  input  from  the  source,  are  recorded 
and  displayed  simultaneously  on  a  square  grid.  The  grid  corresponds  to 
the  locations  of  the  individual  receivers.  For  a  reflecting  surface  parallel 


I 
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to  the  plane  of  the  array,  the  reflection  will  be  observed  as  a  circular 
pattern  centered  on  the  source  location.  The  circles  expand  from  the 
fixed  center  with  increasing  time.  If  the  reflector  is  dipping,  then  the 
center  of  the  pattern  is  shifted  updip  to  a  point  which  is  the  projection 
of  the  image  of  the  source  onto  the  plane  of  ihe  array  (Fig.  3-7). 

The  recorded  signals  from  the  receivers  are  viewed  on  a  screen  as  a 
variable  intensity  display  of  all  the  receiver  outputs  at  a  particular 
time.  As  time  is  varied  by  scanning  through  the  record,  reflections 
(and  other  coherent  events)  will  appear  as  circular  rings  or  arcs  of 
circles,  move  across  die  screen  and  then  disappear. 

The  display  constantly  changes  as  time  is  increased  and  deeper  and  deeper 
reflections  appear.  The  rings  of  reflected  energy  get  thinner  and  move  more 
slowly  as  they  expand  toward  die  edges  of  the  display.  Since  the  rings  of 
reflected  P-wave  energy  will  move  at  a  different  rate  than  those  of  S -waves, 
and  both  will  move  at  different  rates  than  noise,  they  are  distinguishable. 

When  crossed  arrays  are  used  to  simplify  field  procedures,  the  interpretation 

becomes  somewhat  more  complex.  For  a  single  source,  all  rays  spread 

o 

out  radially  from  the  source  to  the  N  receivers.  For  a  crossed  array  of 

2 

N  sources  and  N  receivers,  N  signals  are  still  recorded,  but  the  circular 
symmetry  of  die  single  source  geometry  is  lost.  As  a  result,  the  circular 
patterns  are  distorted  into  ellipses,  whose  eccentricities  depend  on  the 
angles  of  dip  of  die  reflectors. 

Figures  3-8  and  3-9  are  simulated  displays  for  reflectors  dipping  to  the 
left  (Fig.  3-8)  and  to  the  lower  left  (Fig.  3-9).  Individual  frames,  from 
left  to  right,  represent  successive  instants  of  time  at  ten  microsecond 
increments,  with  an  assumed  P-wave  velocity  of  17,  500  ft/sec.  The 
purpose  of  the  simulation  was  to  get  a  preliminary  evaluation  of  the 


Figure  3-8.  Simulated  Displays  for  Reflectors  Dipping  to  the  Left 
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Figure  3-9.  Simulated  Displays  for  Reflectors  Dipping 
to  the  Lower  Left 
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applicability  of  the  3-D  method  to  relatively  small  arrays,  which  might 
be  practical  for  underground  use.  Figures  3-8  and  3-9  simulate  21  x  21 
element  arrays,  in  contrast  to  48  x  48  element  arrays  used  for  seismic 
exploration  by  the  3-D  method.  The  seismic  pulse  was  assumed  to  be  a 
damped  sinusoid.  The  element  spacing  is  one-half  wavelength  at  the 
dominant  frequency  of  the  pulse.  Consequently,  the  simulated  array 
aperture  is  ten  wavelengths  square.  Sixteen  intensity  grey  levels  were 
used. 

Wavefront  distortion  from  causes  listed  above,  which  may  be  expected 
to  severely  degrade  beamforming  techniques  will  merely  cpstort  the 
patterns  observed  with  the  3-D  method.  To  the  extent  that  the  patterns 
are  st:!l  recognizable,  the  reflections  may  be  more  detectable  than  ihey 
would  be  with  the  beamforming  technique.  Also,  faults  on  a  reflection 
will  break  its  circular  pattern,  and  the  location  of  the  break  indicates 
where  the  fault  is.  Again,  using  beamforming,  the  presence  of  a  fault 
will  partly  destroy  the  phase  coherence  of  the  reflected  wavefront,  and 
reduce  the  output  of  the  beamformer  accordingly.  This  represents  an 
advantage  of  the  3-D  method  for  underground  applications,  since  fault 
detection  is  of  high  priority. 
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SECTION  IV 


4.  0  SEISMIC  PROTOTYPE  SYSTEM  DEVELOPMENT 

The  objective  for  the  first  six  months’  phase  of  the  present  program 
has  been  the  design  and  construction  of  a  seismic  system  for  further 
field  experiments  to  evaluate  beamforming  and  3-D  seismic  array 
methods.  The  goal  is  a  prototype  system  for  application  of  the  3-D 
method  which  can  also  provide  recordings  for  subsequent  "playback", 
digitization,  and  digital  processing  in  the  laboratory  for  evaluation  of 
beamforming.  Desirable  characteristics  of  the  system  include: 

(1)  Portability 

(2)  Low  cost 

(3)  Minimum  power  requirements 

(4)  Flexibility 

(5)  Simplicity 

(6)  Maximum  accuracy  and  speed  of  operation 


4. 1  System  Concept 

A  functional  block  diagram  of  the  system  is  shown  in  Figure  4-1.  A 
seismic  recording  system  consisting  of  single  source  transmitter,  single 
receiver  and  optical  recorder  produces  variable  density  film  transparencies 
which  can  be  viewed  in  the  field  using  the  display  system,  or  digitized 
in  the  laboratory. 

The  display  system  contains  a  coherent  fiber  optics  image  converter.  The 
film  format  consists  of  a  variable  density  record  from  each  array  element 
side-by-side  along  one  axis,  with  record  time  along  the  other  axis  (elapsed 


Transmitter  I  Receiver 


Prototype  Seismic /Acoustic  System 


time  relative  to  transmit  pulse).  The  fiber  optics  converter  has  a  row 
of  fiber  bundles  which  samples  all  array  element  records  corresponding  to 
a  particular  time.  The  fiber  bundles  are  rearranged  at  the  other  end  of 
the  device  into  a  square  matrix  corresponding  to  the  seismic  array.  The 
image  observed  at  the  square  end  is  the  array  output  at  a  particular  time 
as  required.  Time  is  scanned  by  moving  the  transparency. 

A  single  record  can  be  displayed  at  a  time  on  the  cathode -ray  tube  (CRT). 
The  Z-axis  (intensity)  of  the  CRT  is  modulated  by  the  receiver  response 
to  a  transmitted  pulse.  A  linear  sweep  signal  deflects  the  CRT  spot  to 
produce  an  intensity -modulated  line  display  of  the  received  signal  versus 
time.  The  line  is  photographed,  producing  a  variable  density  record  of 
one  element  of  the  array.  The  receiver  is  moved  to  the  next  array 
location,  the  film  is  moved  one  record  width,  and  the  process  repeated. 
This  procedure  is  repeated  until  all  array  locations  have  been  recorded. 

Time  synchronization  is  maintained  between  records  by  triggering  the 
sweep  generator  by  the  pulsed  source  transmitter.  Delaying  the  triggered 
sweep  allows  a  selected  time  interval  to  be  recorded. 

The  delayed  sweep  generator  function  was  implemented  in  the  system  by  a 
Biomation  Model  802  Transient  Recorder.  This  commercial  unit  provides 
a  memory  function,  which  allows  nonrepetitive  seismic  sources,  such  as 
hammers  and  other  impact  devices,  to  be  used  for  more  flexibility. 

Polaroid  Type  55P/N  film  was  specified  for  use  in  the  recording  camera, 
as  a  matter  of  convenience  for  field  use.  The  film  provides  both  a 
transparency  and  a  print,  has  sufficiently  fine  grain  in  the  negative  for  the 
required  magnification,  and  good  grey  scale  range. 
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Figure  4-2  shows  the  completed  seismic  recording  system,  less  camera. 
The  configuration  shown  includes  a  borehole  receiver  and  spark  transmitter, 
which  was  developed  as  a  borehole  source.  The  control  unit  on  the  right 
also  has  provision  for  a  piezoelectric  ceramic  (PZT)  transducer  (not  shown) 
which  will  be  evaluated  as  a  surface  source. 

Descriptions  of  construction  details  and  modifications  are  given  in  the 
discussion  to  follow  for  the  optical  recorder,  receiver,  and  display 
subsystem.  The  seismic  source  development  will  be  discussed  in  Section 
5.0. 

4. 2  Optical  Recorder 

The  optical  recorder  consists  of  the  Biomation  Model  802  Transient 
Recorder,  Tektronix  Type  602  Display  Unit,  and  modified  Tektronix  C-27 
Oscilloscopt  Camera  (Fig.  4-3). 

The  802  digitally  stores  1C24  samples  of  the  signal  from  the  receiver 
accelerometer  (also  shown).  The  samples  represent  the  receiver  output 
over  a  time  interval  determined  by  the  sweep  time  and  delay  settings 
chosen.  Thus,  for  example,  it  is  possible  to  gate  out  the  first  millisecond 
of  receiver  output  and  record  the  succeeding  5  milliseconds. 

The  contents  of  the  digital  memory  are  converted  to  a  repetitive  analog 
output  signal  synchronized  with  a  repetitive  two -millisecond  linear  sweep. 
The  Y-output  (signal)  and  X -output  (sweep)  drive  the  Z-axis  (intensity)  and 
Y-axis  (vertical  deflection),  respectively,  of  the  602  CRT  display.  The 
receiver  output  is  feus  made  repetitive  and  displayed  as  an  intensity- 
modulated  line,  just  visible  in  the  center  of  fee  CRT  face  in  Figure  4-3. 
Figure  4-4  shows  a  close-up  of  the  CRT  face. 


A  vertical  slit  aperture,  16  mils  wide,  has  been  added  to  the  CRT  face¬ 
plate.  The  slit  serves  to  provide  a  uniform  record  width.  The  CRT  spot 
size  normally  increases  with  intensity,  which  is  proportional  to  the 
receiver  output  voltage.  To  avoid  variations  in  trace  width  with  signal 
amplitude,  the  slit  was  provided,  together  with  a  circuit  to  provide  a  high 
frequency  triangular  wave  signal  drive  on  the  CRT  X-axis  (Fig.  4-5). 

This  higji  frequency  "dither”  produces  a  broadened  trace  of  uniform 
intensity,  many  times  the  width  of  the  slit.  In  addition  to  providing  a 
uniform  record  width,  this  approach  eliminates  horizontal  drift  as  a 
possible  source  of  error. 

The  CRT  has  only  an  outline  graticule,  so  that  no  lines  will  appear  in  the 
film  record.  (An  external  time -mark  generator  can  be  used,  however,  to 
record  variable -density  time  lines,  when  desired. )  The  outline  graticule 
appears  as  a  fine  line  near  die  top  and  bottom  of  the  slit  which  can  be  used 
as  index  marks  for  initial  vertical  trace  positioning  and  periodic  correction 
of  vertical  drift  with  the  CRT  vertical  position  control,  if  necessary.  The 
phase  coherence  of  array  outputs  can  thus  be  maintained  to  a  high  degree 
of  accuracy. 

4. 2. 1  Camera  Modifications —  The  camera  assembly  consists  of  camera, 
camera  back  and  film  ho.1  Jer,  and  film  advance  mechanism.  The  addition  of 
the  lead  screw  detented  film  advance  was  the  principal  modification,  although 
changes  were  also  made  to  the  camera  back. 

The  recording  film  must  be  precisely  movable  in  increments  slightly  less 
than  one  trace  width  in  order  to  avoid  undesirable  gaps  between  successive 
traces  on  the  film.  This  wais  accomplished  with  a  lead  screw  to  move  the 
camera  back.  The  lead  screw  contaiins  32  threads  per  inch  and  is  detented 
at  quarter-turns  (Fig.  4-6),  so  each  detent  position  represents  a  film  motion 
of  1/128  inches  or  about  0. 0078".  The  camera  magnification  is  1:0.  5,  so 
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that  the  image  of  the  CRT  slit  aperture  is  half  the  width  of  the  slit. 

A  slit  width  of  0.  016  inch  produces  a  0.  0080-inch  record,  which  allows 
0. 0002  inch  record  overlap. 

The  modified  camera  will  record  128  records  per  inch  of  film,  without 
gaps  or  excessive  overlap  between  adjacent  records.  Figure  4-7  is  a 
test  recording,  containing  384  traces.  The  recorded  signal  was  a  sine 
wave,  with  the  frequency  varied  slightly  between  adjacent  traces. 

Resolution  along  the  time  axis  appears  to  be  at  least  100  cycles/sweep. 

For  a  ten  millisecond  sweep,  compatible  with  penetration  up  to  about 
100  feet,  this  resolution  provides  a  frequency  response  of  about  10  KHz. 

A  standard  4x5  Graflok  camera  back  was  modified  with  extender  plates 
to  allow  movement  of  the  film  over  a  total  distance  of  four  inches.  The 
four-inch  allowed  movement  permits  the  recording  of  up  to  512  array 
outputs.  Thus,  maximum  array  dimensions  sure  22  x  22  =  484  elements. 

The  lead  screw  was  made  long  enough  to  accommodate  the  four -inch 
format. 

4. 2. 2  System  Sensitometric  Measurements  —  Densitometer  measurements 
were  made  on  the  developed  film  to  optimize  the  combination  of: 

(1)  Signal  input  levels  to  the  CRT  Z-axis, 

(2)  Camera  aperture,  and 

(3)  Exposure  time. 

The  best  combination  consisted  of : 

(1)  Signsd  level:  0-1.0  volt 

(2)  Aperture:  f/5. 6 

(3)  Exposure  time:  1  second 


A  representative  graph  of  film  density  versus  input  voltage  is  shown  in 
Fig.  4-8.  The  range  of  density  is  such  that  16  gray  levels  were  easily 
distinguishable. 


4.  3  Borehole  Receiver 

The  receiver  was  designed  to  be  used  either  in  a  borehole  or  on  the  surface 
of  the  mine  or  tunnel  opening.  An  accelerometer  is  used  as  the  transducer. 
For  surface  mounting,  the  accelerometer  can  be  cemented  or  used  with  a 
detachable  stud. 

A  pneumatic  method  was  devised  for  clamping  the  accelerometer  to  the 
side  of  a  small  diameter  borehole.  Semi-rigid  plastic  tubing  is  used  to 
position  the  accelerometer  in  the  hole  and  to  inflate  the  rubber  clamping 
hose,  using  a  bicycle  tire  pump. 

The  assembly  is  pictured  in  Fig.  4-9,  inflated  inside  a  1-1/2  inch  diameter 
plastic  tube.  A  replaceable  steel  or  aluminum  button  is  used  to  hold  the 
accelerometer  in  place,  provide  good  mechanical  coupling  to  the  rock 
surface,  and  protect  the  accelerometer.  Insulated  washers  are  used  for 
electrical  isolation. 

Figure  4-10  shows  a  closeup  of  the  device  wit'  the  accelerometer  removed. 
The  rubber  hose  is  clamped  to  a  brass  cylinder,  which  is  hollow  at  one  end 
for  inflating  the  tube. 

4. 3. 1  Accelerometer  --  A  commercially -available  piezoelectric  accelerom¬ 
eter,  Endevco  Model  2220C,  ^as  chosen  as  die  transducer.  Selection 
criteria  included:  1)  small  size,  2)  higi  resonant  frequency,  3)  center-hole 
configuration,  and  4)  sensitivity  (Table  IV- 1). 


Figure  4-8.  Gray  Scale  Measurements  of  Panatomic  and 
Polaroid  P/N  Films 
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Table  IV- 1*  Specifications  for  Model  2220C  Accelerometer 


Charge  Sensitivity 
Voltage  Sensitivity* 
Transducer  Capacitance 


2. 8  pC  /g,  nominal 
3.  0  mV f  g,  nominal 
750  pF,  nominal 


Mounted  Resonant  Frequency  50,  000  Hz,  nominal 


Resistance 


Frequency  Response 
Transverse  Sensitivity 
Amplitude  Linearity 


20,  000  M  Cl ,  minimum  at  +72°  F 
5000  MCI,  minimum  at  +350°  F 

2  Hz  to  10,  000  Hz 

5%  maximum,  2%  special  selection 

Sensitivity  increases  approximately  l%per 
500g,  0  to  5000g 


with  200  pF  external  capacitance 


V**r  Vt-v-h; 


0 
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) 

The  high  frequency  response  of  a  piezoelectric  accelerometer  is  determined  ) 

by  its  mechanical  resonance.  It  can  be  represented  as  an  undamped 

resonant  spring-mass  system,  and  provides  a  flat  response  to  acceleration  \ 

inputs  up  to  about  one -fifth  its  resonant  frequency,  with  minimum  phase 
shift.  For  pulse  measurements,  accelerometer  overshoot  and  ringing  are  j 

minimized  if  the  natural  period  is  less  than  one -fifth  the  pulse  duration,  or 

less  than  about  one-third  Ihe  pulse  rise  time.  Since  we  are  concerned  with  ) 

) 

pulse  durations  of  about  100  microseconds,  and  frequencies  below  10  KHz, 

the  50  KHz  resonant  frequency  of  the  2220C  accelerometer  was  considered  i 

adequate.  J 

The  low  frequency  response  of  a  piezoelectric  accelerometer  is  determined  ] 

when  used  with  voltage  amplification  by  the  RC  time  constant  of  the  total 

capacitance  (transducer,  cable  and  input  capacitance  of  the  amplifier)  along  !  j 

with  the  input  resistance  of  the  amplifier.  When  charge  amplification  is 

used,  the  low  frequency  response  is  determined  primarily  by  the  amplifier  |  I 

'  } 

response. 

4.  3. 2  Amplifier  —  A  bandpass  amplifier  was  incorporated  into  the  receiver,  J 

whose  primary  purpose  is  to  provide  high  frequency  roll-off,  so  that 

accelerometer  resonance  response  is  reduced,  and  to  provide  a  low  frequency  J 

response  which  will  be  independent  of  cable  capacitance,  and  therefore  the 

length  of  cable  used.  The  charge  amplifier  circuit,  with  low  and  high  ) J 

frequency  cutoff  points  at  2  KHz  and  8  KHz  and  6  db/octave  rolloff,  is  shown 

schematicaUy  in  Fig.  4-11.  j 

] 

] 

i  | 

'I 

I 
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SECTION  V 


5.  0  SEISMIC  SOURCE  EVALUATION  AND  DEVELOPMENT 

A  substantial  effort  has  been  directed  in  this  phase  toward  the  evaluation 
of  possible  seismic  sources  for  use  on  the  rock  surface,  as  well  as  in  a 
borehole.  A  series  of  measurements  were  taken  by  the  University  of 
Minnesota  under  subcontract  for  a  variety  of  surface  sources,  primarily 
of  impact  type.  A  summary  of  these  measurements  is  contained  below. 
In  addition,  a  pulsed  piezoelectric  source,  developed  by  the  U.  S.  Army 
for  tunnel  detection  and  recently  tested  as  a  seismic  source  in  a  coal 
mine,  was  acquired  for  evaluation  as  source  transmitter  in  planned 
underground  data  acquisition.  Finally,  an  underwater  spark  transmitter 
was  developed  to  serve  as  a  downhole  source  for  array  experiments. 


5. 1  Summary  of  Impact  Source  Measurements 

A  series  of  measurements  was  taken  to  obtain  characteristics  of  seismic 
waveforms  produced  by  a  variety  of  impact-type  sources  including: 

(1)  A  hammer,  directly  onto  the  rock  surface, 

(2)  Hammer,  onto  pins  embedded  in  the  rock, 

(3)  Hammer,  onto  rods  in  contact  with  the  rock,  both  on  the 
surface  and  at  the  bottom  of  a  hole,  and 

(4)  Explosive-  and  gas -driven  projectiles. 
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Waveforms  were  detected  by  piezoelectric  accelerometers  and  converted 
to  permanent  stripchart  records  {with  special  instrumentation  to  provide 
frequency  response  to  100  KHz)  and  oscilloscope  photographs.  The 
waveform  was  passed  through  a  filter  which  permitted  variable  bandpass 
filtering.  The  accelerometer  and  source  placements  were  selected  to 
yield  results  for: 

a)  Transmission  through  a  concrete  block, 

b)  Propagation  along  the  surface, 

c)  Detection  within  a  drill  hole. 

Accelerometer  tests  within  the  drill  hole  were  carried  out  using  two  methods 
of  holding  the  accelerometer  against  the  side  of  the  hole:  mechanical  and 
pneumatic.  The  results  indicated  no  difference  in  rise  time,  amplitude,  or 
total  waveform. 

Frequency  content  of  the  waveforms  did  not  change  drastically  among  the 
various  impact -type  sources.  Meat  of  the  energy  is  in  the  range  5-20  KHz, 
with  substantial  contributions  up  to  30-40  KHz  and  negligible  energy  below 
2  KHz.  Impacts  onto  rods,  and  particularly  impacts  onto  embedded  pins,  do 
give  somewhat  more  high -frequency  components  than  direct  impact  onto  the 
rock. 

The  repeatability  of  these  waveforms  is  reasonably  good  for  frequencies  of 
less  than  15-20  KHz,  but  poor  for  higher  frequencies. 

One  measure  of  the  frequency  characteristics  of  the  source  waveforms  is 
the  rise  time  (zero  to  peak)  of  the  direct  P  wave.  This  measure  has  the 
advantage  that  it  will  be  unaffected  by  reflections  and  other  disturbances. 

It  was  found  that: 


(1)  Hammer  impact  directly  onto  the  concrete  surface  gives  a 
rise  time  of  the  order  of  30-40  microseconds. 

(2)  Impact  onto  an  embedded  pin  reduces  the  rise  time  to  10-15 
microseconds. 

(3)  Impacts  onto  a  small  center  punch  and  a  nailset  gives  rise 
times  of  about  25  and  15  microseconds,  respectively. 

(4)  Impacts  onto  a  star  drill  and  a  large  center  punch  gives 
rise  times  of  about  35  and  30  microseconds,  respectively. 

(5)  Impacts  onto  the  steel  and  aluminum  rods  gives  rise  times 
of  about  20  and  15  microseconds,  respectively. 

The  explosive  stud  driver  with  projectile  yields  two  distinct  arrivals. 
Arrival  time  for  the  second  varies  with  explosive  charge  size  but  fails  in 
the  range  350-500  microseconds.  Amplitudes  are  5-10  times  as  great  as 
for  the  first  arrival.  A  consideration  of  possible  explanations  leads  to 
the  conclusion  that  the  smaller  first  arrival  is  from  the  explosive  gases 
and  the  second  from  the  projectile  impact. 

It  was  concluded  that  impact  sources  have  the  following  advantages  as 
possible  seismic  sources  for  the  excavation  seismology  program: 

(1)  They  are  rapid,  cheap,  light,  flexible. 

(2)  They  bypass  problems  associated  with  rough  or  fractured 
rock  surfaces. 


(3)  They  are  rich  in  frequencies  at  the  range  of  interest  for 
this  program,  say  5-20  KHz. 

(4)  Rise  times  are  quite  satisfactory,  especially  for  impacts 
onto  rods  or  embedded  pins. 

(5)  Energy  output  is  good. 

(6)  Waveform  repeatability  is  adequate  in  the  frequency  range 
of  interest,  although  poor  at  higher  frequencies  (which 
could  be  filtered  out  if  desired). 

(7)  Source  location  can  be  precisely  defined  for  array  studies, 
especially  for  impacts  onto  rods  or  embedded  pins. 


5.1  2 3  Underwater  Spark  Source  Development 

Requirements  for  a  high  energy  seismic  impulse  source,  which  could  be 
used  at  depths  up  to  15-20  feet  in  a  small  diameter  borehole,  led  to  the 
development  of  an  underwater  spark  source.  The  source  was  required  to 
be  repetitive,  portable,  and  have  a  pulse  duration  of  not  more  than  about 
100  microseconds,  with  substantial  energy  output  at  frequencies  below  10  KHz 

Other  possibilities  which  were  considered,  but  appeared  less  advantageous 
than  an  underwater  spark  included: 

(1)  Detonators  (blasting  caps), 

(2)  Air  gun, 

(3)  Piezoelectric  or  magnetos  trie  tive  transducer. 


H 

0 

E 
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(4)  Propane/ air  exploder,  and 

(5)  Electromechanical  or  mechanical  impactor. 

Blasting  caps,  although  fairly  inexpensive  and  convenient,  would  release 
gaseous  explosion  products  into  the  water-filled  hole.  Since  we  would  want 
to  fire  repeatedly  several  h-imdred  times  in  the  same  hole,  the  pulse 
characteristics  likely  would  change.  Also,  firing  time  delay  uncertainties 
would  adversely  affect  array  output  time  alignment. 

Air  guns,  which  have  been  used  as  sources  for  offshore  seismic  exploration, 
operate  on  the  sudden  release  of  a  small  volume  of  high  pressure  air.  The 
output  pulse  of  existing  air  guns  was  too  long  in  duration  and  low  in  frequency 
for  the  underground  seismic  problem,  and  it  did  not  appear  feasible  to 
change  it  substantially.  Again,  time  synchronization  would  also  be  a  problem. 

Piezoelectric  or  magnetostrictive  transducers  have  excellent  repeatability 
and  can  be  precisely  triggered.  However,  size  was  considered  to  be  a 
problem.  The  pulse  energy  produced  is  directly  related  to  the  size  of  the 
device.  For  a  small  diameter  hole,  the  allowable  size  would  not  produce 
sufficient  energy  output. 

A  propane/air  exploder  was  a  promising  alternative.  A  very  high  output 
energy  pulse  of  short  duration  could  be  achieved,  and  the  hole  would  not  need 
to  be  filled  with  water  for  coupling  the  seismic  pulse  to  the  rock.  A 
preliminary  design  of  such  a  device  was  obtained;  however,  the  propane 
source  would  have  to  be  retrieved  from  the  hole  after  each  shot  for  reloading 
and  replacement  of  a  frangible  disc,  Since  this  operation  would  be  undesirable, 
no  further  development  of  this  source  was  attempted. 


'  **  v 
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Mechanical  impactor  concepts  were  severely  limited  by  the  limited  hole 
size,  and  again,  time  synchronization  would  be  a  problem  with  mechanical 
devices. 

5.  2. 1  The  Electrohydraulic  Effect  --An  underwater  spark  produces  a 
pressure  pulse  which  has  found  application  as  a  sonic  source  (Ref.  1) 
and  has  been  investigated  for  rock  fracturing  and  comminution  (Ref.  4). 

The  underwater  spark  is  produced  by  the  sudden  discharge  of  capacitively 
stored  electric  energy  across  a  submerged  electrode  gap.  Initiation  of  the 
spark  creates  a  plasma  and  vapor  channel  between  the  two  electrodes. 
Extremely  high  pressures  are  produced  within  the  channel,  which  are 
transmitted  to  the  surrounding  water.  The  channel  then  expands  as  a 
bubble,  compressing  the  water  until  the  vapor  pressure  in  the  bubble  falls 
below  the  hydrostatic  pressure.  The  pressure  pulse  transmitted  to  the 
water  is  similar  to  that  from  an  explosion,  rising  rapidly  and  decaying  more 
slowly. 


5. 2. 2  Bubble  Pulsations  —  The  water  in  the  immediate  region  of  the 
expanding  bubble  has  a  large  outward  velocity  and  the  diameter  of  the  bubble 
increases  rapidly.  As  the  pressure  falls  below  the  equilibrium  level,  the 
motion  is  reversed  abruptly,  and  the  bubble  contracts. 

As  the  bubble  expands  and  contracts,  a  pressure  wave  propagates  radially 
outward.  From  incompressible  fluid  theory,  it  is  known  that  pressure 
depends  on  the  square  of  the  rate  of  bubble  expansion  and  contraction.  This 
rate  is  greatest  when  the  bubble  is  at  the  point  of  minimum  volume.  Thus, 
pressures  are  significant  only  near  the  point  of  maximum  contraction.  The 
peak  pressure  in  the  bubble  pulse  is  no  more  than  ten  to  twenty  percent  of 
that  in  the  main  shock  wave  but  the  duration  is  much  greater  (Ref.  2). 
Usually,  only  the  first  bubble  pulse  is  of  practical  importance. 
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5. 2. 3  Circuit  Analysis  —  Figure  5-1  contains  a  functional  block  diagram 
of  the  spark  source.  The  source  is  basically  a  series  network  of  a 
charged  capacitor,  fast  discharge  switch,  connecting  wires,  and  the 
electrode  gap.  Figure  5-2  is  the  equivalent  electrical  circuit  where: 


Voltage  across  charged  capacitor 
Capacitance  of  storage  capacitor 
Inductance  of  storage  capacitor 
Inductance  of  switch  capacitor 
Resistance  of  switch  capacitor 
Capacitance  of  cable  capacitor 
Inductance  of  cable  capacitor 
Capacitance  of  electrode  gap 
Inductance  of  electrode  gap 
Resistance  of  electrode  gap 


The  shunt  capacitances  associated  with  the  cable  and  electrode  gap  are 
negligible,  resulting  in  the  simplified  circuit  of  Fig.  5-3.  Here: 


L  =  Total  inductance  =  L„  +  L„  +  L# 

C  8  v 

R  =  Total  resistance,  not  including  gap  =  Rg  + 

For  a  sudden  closure  of  the  switch,  the  equation  describing  the  discharge 
current  is: 

L  &  +  RI  +  5  J  I  dt  =  o. 

whose  solution  (underdamped  case)  for  I  (t  =  0)  =  0  is 
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I(t)  =  -rn~  e'bt  sin  u>  t 
w„i-<  o 


where 


“o  =  y 


1 

LC 


4L 


tj—  =  frequency  of  damped  sinusoid 


R 


b  =  —  =  damping  factor . 


The  maximum  value  of  current  is  determined  from  the  condition 


sL 

dt 


=  0. 


I  =  I. 


max 


Thus, 


max  a)  L 
o 


e  ‘b‘max  sin  a  t 

o  max 


max 


1 


tan 


-1 


u> 


=  time  of  peak  current. 


5. 2. 4  Energy  Considerations  —  The  total  energy  stored  in  the  capacitor 
is  equal  to  =  y  CVQ^.  The  energy  dissipated  by  the  spark  is  given 

by 


E 


jjf 

spark  "  R  +  Rg  Etotal 
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Maximum  efficiency  requires  that  R  <<  R  .  The  spark  resistance  R 

w 

is  determined  by  the  maximum  electrode  gap  which  will  produce  breakdown 
at  the  voltage  V  . 

5. 2.  5  Electrode  Design  —  The  coaxial  electrode  design  which  evolved 
diming  the  spark  source  development  is  shown  in  Fig.  5-4.  The  high 
voltage  electrode  is  the  center  screw.  The  electrode  exposed  area  is 
minimized  to  eliminate  leakage  paths  as  much  as  possible.  The  electrode 
gap  is  filled  with  a  polyester-glass  washer  to  provide  mechanical  rigidity 
and  eliminate  fouling.  Sparking  thus  occurs  around  the  periphery  of  the 
electrode  gap. 

The  maximum  gap  size  for  voltages  up  to  2.  5  Kv  in  tap  water  is  about  0. 025 
inch.  We  found  that  rubber  "O'1  rings,  placed  over  fee  gap  and  held  in 
place  by  fee  linen  bakelite  cap,  would  effectively  suppress  the  formation  of 
a  bubble  pulse.  Figure  5-5  shows  the  electrode  housing,  partially 
disassembled.  The  electrode  structure  can  be  simply  unscrewed  for  replace¬ 
ment. 


This  design  was  found  to  be  highly  rugged  and  to  provide  reliable,  repeatable 
sparking. 


r*  rr*  *"****Ax*«1U? 
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The  frequency  is  estimated  from  the  time  difference  t  between  successive 
current  peaks.  The  damping  factor  is  determined  from  the  ratio  of  current 
values: 


Ij  g-bt  max 

In  c  In  _b  (t  +T) 
I  e  max 

c* 


-  bt 


or. 


b 


R 

2L 


Also, 


1 

LC 


+  b 


+ 


b 


2 


Figure  5-6  is  the  current  waveform  with  (a)  the  electrode  gap  shorted,  to 
measure  circuit  resistance,  and  <b)  during  sparking  with  the  electrode  gap 
immersed  in  tap  water.  The  value  of  t,  from  Figure  5-6  (a>  is  approximately 
30  microseconds.  Calculated  values  of  circuit  inductance  and  resistances, 
with  shorted  gap  and  during  sparking  are  summarized  in  Table  V-l. 

The  experimental  results  indicate  that  approximately  half  of  the  total  resis¬ 
tance  is  due  to  the  gap;  therefore,  half  the  energy  is  dissipated  in  the  spark. 
Of  this,  perhaps  half  is  radiated  as  acoustic  energy  in  the  shock  wave.  The 
total  energy  for  C  =  15  microfarads,  VQ  =  2.  5  kilovolts  is  about  48  joules. 

The  estimated  useful  seismic  energy  per  pulse,  therefore,  is  12  joules 
(1  joule  =  1  watt-second). 


Figure  5-6.  Current  Waveforms  for: 

(a)  electrode  gap  shorted 

(b)  during  sparking 

(One  division  =  2  microseconds) 
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5.  2.  7  Energy  Spectrum  of  Pressure  Pulse  —  The  pressure  pulse  from  an 
underwater  spark  is  approximately  a  decaying  exponential: 

r  -mt 

I-  =  P  e 
o 

where  the  decay  constant  m  is  approximately  the  same  as  the  discharge  cir¬ 
cuit  time  constant,  to  within  10%  (Ref.  1). 

The  acoustic  energy  contained  in  the  pressure  pulse  is 

EP  ■  p(t)l2  dt 


where  pc  is  the  acoustic  impedance  of  wate 
spectrum  is  given  by  the  Fourier  Theorem  at 


2P 


Ep(l)  = 


DC 


*  *• 


i 

5  9  J 

rn  +  (2ttfr 


e  corresponding  frequency 


The  output  energy  for  the  spark  source  is  seen  to  decrease  with  frequency. 
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At  f  =  10  kHz,  assuming  msb  =  bxlO  sec  ,  the  output  spectrum  is  down 
by  a  factor  of  two. 


5.  2.  8  Spark  Source  Circuitry  —  Figure  5-7  is  a  schematic  diagram  of  the 
high  voltage  supply,  storage  capacitor  and  switch.  Triggering  of  the  pulse 
can  be  done  either  manually  by  pushbutton,  or  remotely.  An  autotransiormer 
and  meter  have  been  incorporated  to  allow  the  discharge  voltage,  V  ,  to  be 
continuously  adjustable.  A  synchronizing  pulse  is  also  provided  to  trigger  the 
Biomation  recorder.  The  circuit  is  enclosed  in  a  portable  case,  together 
with  triggering  circuitry  for  a  piezoelectric  transducer,  which  will  be  de¬ 
scribed  in  Section  5.  3.  The  case  and  front  panel  are  shown  in  Figure  4-2. 

In  Figure  5-8,  the  panel  has  been  removed. 


Figure  5-8.  Control  Module  with  Front  Panel  Lowered 


The  discharge  circuit  includes  about  15  feet  of  RG  8/U  cable,  connecting  the 
spark  gap  with  the  charging  circuit.  The  total  resistance  of  the  discharge 
circuit,  including  the  cable,  was  derived  in  5ection  5.2.  6  as  0.  09  ohms.  The 
charging  time  is  approximately  five  seconds. 


Table  V-l.  Measured  Circuit  Parameters 


Gap  Shorted 

Underwater  Spark 

t,  sec 

30  x  10‘6 

30  x  10'6 

!1 

2.  5 

6.0 

h 

6 

0.  9 

1  3 

,  -1 
b,  sec 

3.0  x  104 

6.0  x  104 

L,  microhenries 

1.5 

1.  5 

Resistance,  ohms 

.  R  =  0.09 

R  +  Rg  =  0.  18 

5.  3  Piezoelectric  Transducer 

A  piezoelectric  transducer,  similar  to  one  developed  by  the  U.  S.  Army  for 
seismic  tunnel  detection  (Ref.  6)  and  previously  tested  by  the  U.S.  Bureau  of 
Mines  in  underground  coal  mines,  is  an  attractive  possibility  as  a  seismic 
surface  source  in  hard  rocks. 

The  transducer,  pictured  in  Figure  5-9,  is  based  upon  the  rapid  discharge  of 
an  electrical  charge  stored  in  the  capacitive  ceramic  elements  of  the  transducer. 
The  transducer  is  charged  slowly  from  a  high  voltage  source,  and  the  impulsive 
discharge  provides  a  high-peak-power  mechanical  output  pulse. 

The  capacitance  of  the  transducer  is  0.  30  microfarads  and  the  operating  vol¬ 
tage  is  3  kilovolts  for  a  discharge  energy  per  pulse  of  1,  35  joules.  A  series 
inductance  of  6  millihenries  is  included  to  provide  a  pulse  rise  time  of  47.  5 
microseconds,  and  a  series  resistor  provides  optimum  damping. 


5-20 


For  this  program,  we  provided  a  triggering  circuit  and  DC  power  supply  to 
operate  the  piezoelectric  transducer.  The  circuit  schematics  were  included  in 
Figure  5-7.  The  electronics  are  housed  with  the  spark  source  electronics,  so 
that  either  source  may  be  used  interchangeably. 
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